Evaporation of 2-Dimensional Black Holes 
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We present a detailed analysis of results from a new study of the quantum evaporation of Callan-Giddings- 
Harvey-Strominger (CGHS) black holes within the mean-field approximation. This semi-classical theory in- 
corporates back reaction. Our analytical and numerical calculations show that, while some of the assumptions 
underlying the standard evaporation paradigm are borne out, several are not. One of the anticipated properties 
we confirm is that the semi-classical space-time is asymptotically flat at right future null infinity, X^, yet incom- 
plete in the sense that null observers reach a future Cauchy horizon in finite affine time. Unexpected behavior 
includes that the Bondi mass traditionally used in the literature can become negative even when the area of the 
horizon is macroscopic; an improved Bondi mass remains positive until the end of semi-classical evaporation, 
yet the final value can be arbitrarily large relative to the Planck mass; and the flux of the quantum radiation at 
is non-thermal even when the horizon area is large compared to the Planck scale. Furthermore, if the black 
hole is initially macroscopic, the evaporation process exhibits remarkable universal properties. Although the 
literature on CGHS black holes is quite rich, these features had escaped previous analyses, in part because of 
lack of required numerical precision, and in part due to misinterpretation of certain properties and symmetries 
of the model. Finally, our results provide support for the full quantum scenario recently developed by Ashtekar, 
Taveras and Varadarajan, and also offer a number of interesting problems to the mathematical relativity and 
geometric analysis communities. 
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I. INTRODUCTION 

Although literature on the quantum nature of black holes 
is vast, many important questions on the dynamics of their 
evaporation still remain unanswered. This is true even for 2- 
dimensional models, introduced some twenty years ago [1-3]. 
In this paper, as a follow up to [4], we present a detailed anal- 
ysis of the semi-classical dynamics of 2-dimensional, Callan- 
Giddings-Harvey-Strominger (CGHS) black holes [1], using a 
combination of analytical and high precision numerical meth- 
ods. This model has been studied extensively in the past [2]. 
Yet we find new and rather remarkable behavior for classes of 
black holes where the collapse is prompt and the Arnowitt- 
Deser-Misner (ADM) mass is large: As these black holes 
evaporate, after an initial transient phase, dynamics of vari- 
ous physically interesting quantities flow to universal curves. 
In addition, this analysis brings out certain unforeseen limita- 
tions of the standard paradigm that has been used to discuss 
black hole evaporation for some two decades. The overall sit- 
uation bares interesting parallels to the discovery of critical 
phenomena at the threshold of gravitational collapse in clas- 
sical general relativity [5, 6]; at the time it was assumed that 
the spherically symmetric gravitational collapse problem was 
essentially "solved", yet all earlier work had missed this very 
rich and profound property of gravitational collapse. Also, it 
served as a new example of rather simple, universal behavior 
that dynamically emerges from a complicated system of non- 
linear partial differential equations. Here, though the physical 
scenario is different, a similar unexpected universal behavior 



arises. This universality strongly suggests that, although the 
S-matrix is very likely to be unitary, all the information in the 
infalling matter will not be imprinted in the outgoing quan- 
tum state. As discussed in detail in the paper, this mismatch 
between unitarity and information recovery is a peculiarity of 
2 dimensions. 

Our investigation is carried out within the mean-field ap- 
proximation of [7] in which the black hole formation and 
evaporation is described entirely in terms of non-linear, partial 
differential equations (PDEs). Thus there will be no Hilbert 
spaces, operators or path integrals in this paper. ' The focus 
is rather on the consequences of the geometrical PDEs gov- 
erning black hole evaporation in the semi-classical regime, 
and the intended audience is mathematical and numerical rel- 
ativists. Therefore we will start with an introductory review 
of the CGHS model — the classical field theory in Sec. I A and 
semi-classical corrections in Sec. IB — and then summarize 
the results and outline the rest of the paper in Sec. I C. 



A. The Classical CGHS Model 

Consider first the spherical collapse of a massless scalar 
field / in 4 space-time dimensions. Mathematically, it is con- 
venient to write the coordinate r which measures the phys- 



ical radius of metric 2-spheres as r 



/k where k is a 



constant with dimensions of inverse length. The space-time 
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' This paper complements a companion article [8] which is primarily de- 
voted to quantum issues, particularly that of potential information loss, and 
another companion article [9] devoted to the details of numerical simula- 
tions. 
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metric '^gab can then be expressed as 

Sab = gab+r^Sab ■■= gab + —Sab, (1-1) 

where Sab is the unit 2-sphere metric and g^^ is the 2-metric 
in the r-t plane. In terms of these fields, the action for this 
Einstein-Klein-Gordon sector can be written as 

S{g, ^, /)-^ ^Jd'x^\e~''f' (R + 2V>Va(/. 

StT G4 

+ 2e-^^K^) - i/d2.T/y e-^WVa/ (1.2) 

where G4 is the 4-dimensional Newton's constant, V is the 
derivative operator and R the scalar curvature of the 2-metric 
gal,- The CGHS model, by contrast refers to gravitational col- 
lapse of a scalar field in 2 space-time dimensions. The grav- 
itational field is now coded in a 2-metric g^,^ and a dilaton 
field (j), and the theory has a 2-dimensional gravitational con- 
stant G of dimensions [ML]~^ in addition to the constant k 
of dimensions [L]~^ (k^ is sometimes regarded as the cosmo- 
logical constant).- The CGHS action is given by [1]: 

+4^)-^Jd\^\WyWa.f. (1.3) 

Note that the two actions are closely related. The only differ- 
ence is in some coefficients which appear bold faced in (1.2). 
This is why one expects that analysis of the CGHS model 
should provide useful intuition for evaporation of spherically 
symmetric black holes in 4 dimensions. 

On the other hand, the two theories do differ in some impor- 
tant ways and we anticipate that only certain aspects of uni- 
versality will carry over to 4-dimensions. These differences 
are discussed in section IV. Here, we only note one: since the 
dilaton field does not appear in the scalar field action of (1.3), 
dynamics of / decouples from that of the dilaton. Now, since 
our space-time is topologically M^, the physical 2-metric g^j, 
is conformally flat. We can thus fix a fiducial flat 2-metric 77"'' 
and write g"'' = fJry"'', thereby encoding the physical geom- 
etry in the conformal factor il and the dilaton field (j). Next, 
since the wave equation is conformally invariant, 

□ (g)/ = ^ □(,)./ = 0- (1-4) 

/ is only subject to the wave equation in the fiducial flat space 
which can be easily solved, without any knowledge of the 
physical geometry governed by (fJ, 0). This is a key simplifi- 
cation which is not shared by the scalar field f in the spheri- 
cally symmetric gravitational collapse described by ( 1.2). De- 
note by the advanced and retarded null coordinates of 77 so 
that rjab — 2c)(a2^ db)Z^ . Then a general solution to (1.4) on 
the fiducial Minkowski space rj) is simply 

+ (1.5) 



^ In this paper we set c = 1 but keep Newton's constant G and Planck's 
constant h free. Note that since Gh is a Planck number in 2 dimensions, 
setting both of them to 1 is a physical restriction. 



singularity y' 




Figure 1. Penrose diagram of the CGHS black hole formed by the 
gravitational collapse of a left moving field /+ . The physical space- 
time is that part of the fiducial Minkowski space which is to the past 
of the space-like singularity. 



where f± are arbitrary well behaved functions of their argu- 
ments. In the classical CGHS theory, one sets /_ = and 
focuses on the gravitational collapse of the left moving mode 
/-).. As one might expect, the true degree of freedom lies only 
in /+, i.e., fj^ completely determines the geometry. But in 
the classical CGHS model, there is a further unexpected sim- 
plification: the full solution can be expressed as an explicit 
integral involving /+! 

For later purposes, following [7], let us set 

$ := e-^"^ 

and introduce a new field via 

6 = 0"^$ so that g''^ ^ Q^^ ^ri"^ 

Then the geometry is completely determined by the pair of 
fields 0, $, and field equations obtained by varying (1.3) can 
be solved to express 8, $ directly in terms of /+. The result- 
ing expressions for 6 and $ are simpler in terms of 'Kruskal- 
like' coordinates given by 

= e''''^ , and k,x~ = -e~'^''~ . (1.6) 

Given any regular /+, the full solution to the classical CGHS 
equations can now be written as 

O = —K^X'^ X^ 

* = e-|/;"dx+/f {du/dx+r-- (1.7) 

Note that, given any regular /, the fields (0, $) of (1.7) that 
determine geometry are also regular everywhere on the fidu- 
cial Minkowski manifold M°. 

How can the solution then represent a black hole ? It turns 
out that, for any regular /+, the field $ of ( 1 .7) vanishes along 
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a space-like line Along tg then, vanishes, whence 
the covariant metric g^^;, fails to be well-defined. It is easy 
to verify that the Ricci scalar of g^^ diverges there. This is the 
singularity of the physical metric g. The physical space-time 
(M, g_^i^) occupies only that portion of M° which is to the past 
of this singularity (see Fig. 1). 

But does £s represent a black hole singularity? It is easy to 
check that (A/, g^^) admits a smooth null infinity I which has 
4 components (because we are in 2 space-time dimensions): 
Xj^ and coincide with the corresponding 1^ ^nd I^^ of 
Minkowski space-time (M°, ry) while and 1^ are proper 
subsets of the Minkowskian II"*" andZ^^. Nonetheless,!]^ is 
complete with respect to the physical metric g^^^ and its past 
does not cover all ofM. Thus, there is indeed an event horizon 
with respect to hiding a black hole singularity. However, 
unfortunately is not complete with respect to g^f,. There- 
fore, strictly speaking we cannot even ask^ if there is an event 
horizon — and hence a black hole — with respect to ! For- 
tunately, it turns out that for the analysis of black hole evap- 
oration — and indeed for the issue of information loss in full 
quantum theory — only is relevant. To summarize then, 
even though our fundamental mathematical fields (6,$) are 
everywhere regular on full M°, a black hole emerges because 
physics is determined by the Lorentzian geometry of g. 

Although a black hole does result from gravitational col- 
lapse in the CGHS model, it follows from the explicit solu- 
tion (1.7) that one does not encounter all the rich behavior 
associated with spherical collapse in 4 dimensions. In partic- 
ular there are no critical phenomena [5, 6], essentially because 
there is no threshold of black hole formation: a black hole re- 
sults no matter how weak the infalling pulse /+ is. However, 
as remarked in the beginning of this section, the situation be- 
comes more interesting even in this simple model once one al- 
lows for quantum evaporation and takes into account its back 
reaction. 



B. The Semi-Classical CGHS Model 

To incorporate back reaction, one can use semi-classical 
gravity where matter fields are allowed to be quantum but 
geometry is kept classical. In this paper we will implement 
this idea using the mean field approximation of [7, 8] where 
one ignores the quantum fluctuations of geometry — i.e., of 
quantum fields (G), $) — but keeps track of the quantum fluc- 
tuations of matter fields. The validity of this approximation 
requires a large number of matter fields with i = 1, . . . N 
(whence it is essentially the large N approximation [1, 2]). 
Then, there is a large domain in space-time where quantum 
fluctuations of matter can dominate over those of geome- 
try. Back reaction of the quantum radiation modifies classi- 
cal equations with terms proportional to NGh. However, dy- 
namics of the physical metric g is again governed by PDEs 



singularity y' '-.last ray 




Figure 2. Penrose diagram of an evaporating CGHS black hole in the 
mean field approximation. Because of quantum radiation the singu- 
larity now ends in the space-time interior and does not reach or 
(compare with Fig. 1.) Space-time admits a generalized dynami- 
cal horizon whose area steadily decreases. It meets the singularity at 
its (right) end point. The physical space-time in this approximation 
excludes a future portion of the fiducial Minkowski space (bounded 
by the singularity, the last ray and the future part of the collapsing 
matter). 



on classical fields, {Q, $), which we write without an under- 
bar to differentiate them from solutions (8, $) to the classical 
equations (Nh ~ 0). In the domain of applicability of the 
mean-field approximation, they are given by expectation val- 
ues of the quantum operator fields: 9 = (O) and $ = ($). 
The difference from the classical case is that the coefficients 
of the PDEs and components of the metric gab now contain h. 

For any given finite N, there is nonetheless a region in 
which the quantum fluctuations of geometry are simply too 
large for the mean field approximation to hold. This is re- 
flected in the fact that a singularity persists in this approxima- 
tion, although it is now weakened. The field $ now assumes 
a non-zero value NGh/12 at the singularity whence g"'^ is in- 
vertible there [2]. Furthermore it is C° across this singularity 
but not . Finally, because of back-reaction, the strength of 
the singularity diminishes as the black hole evaporates and the 
singularity ends in the interior of space-time; in contrast to the 
classical singularity, it does not reach (see Fig. 2). It is the 
dynamics of gat that exhibit novel, universal features. 

It turns out that the fundamental equations of the mean-field 
theory (and in fact also of the full quantum theory [7, 8]) ad- 
mit a scaling symmetry, which was discovered independently 
by Ori [11]. This has important consequences, which to our 
knowledge have not been fully appreciated before."* It natu- 
rally leads to certain physical quantities that have universal 



^ Even in 4-dimensions, the black hole region is defined as B := M \ 
J~ (I^) provided I'^ is complete. If we drop the completeness require- 
ment, even Minkowski space would admit a black hole! See, e.g., [10]. 



In [12] for example it was noted that A'^ could be "scaled out" of the prob- 
lem and that the results are "qualitatively independent of N". However, as 
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properties during and at the end of the quantum evaporation 
of large black holes. Unraveling this unforeseen behavior 
requires the use of an appropriate analytical formulation of 
the problem, and high precision numerical solutions (both in 
terms of requiring small truncation error and using the full 
range of 16-digit double precision floating point arithmetic) 
— for details, see the companion papers [8, 9]. These uni- 
versal properties have their origin in the PDEs governing the 
dynamics. Therefore, our numerical results lead to a series 
of interesting conjectures in mathematical relativity and ge- 
ometric analysis (discussed in section IV). In addition our 
detailed analysis shows, quite surprisingly, that several of the 
standard assumptions that have driven the quantum evapora- 
tion paradigm are incorrect. As a consequence, results of this 
paper also play a significant role in resolving the information 
loss puzzle in the quantum theory [7, 8]. 



C. Summary of New Results 

To summarize the new results, let us first recall the standard 
paradigm. Literature on the quantum evaporation of CGHS 
black holes uses a certain definition of Bondi mass M'^'i.fA:. 

rsonai 

Essentially every preceding paper assumed that: i) The semi- 
classical approximation is excellent until the horizon shrinks 
to Planck size; ii) Throughout this long phase, A/g^fj^j is non- 
negative and the process is quasi-static; iii) Consequently, dur- 
ing this phase the quantum flux at is given by the Hawking 
thermal flux and the semi-classical approximation holds; and 
iv) At the end of this phase the Bondi mass is also of Planck 
size. It is then difficult to imagine how purity of the incoming 
quantum state could be preserved in the outgoing state. How- 
ever, our results show that several features of this scenario fail 
to be borne out by detailed calculations in the semi-classical 
theory. In particular, in section III we will show the following 
results for a prompt collapse of data with large ADM mass: 

• The traditional Bondi mass, Mul.f-i:, in fact becomes 
negative (and large) even while the horizon area is 
macroscopic. 

• The definition of Ai^^^j is taken directly from the 
classical theory where the black hole is static. Now, 
in 4-dimensions one knows [13] that the formula for 
the Bondi mass has to be modified in non-stationary 
space-times (from / v]/o d^V^ to /(^g - o-^) d'^V). In- 
deed if one were to ignore this modification, one would 
find that neither the Bondi mass nor the Bondi flux 
is always positive. By mimicking the 4-d procedure, 
Ashtekar, Taveras and Varadarajan (ATV) [7] had pro- 
posed a quantum corrected Bondi mass, A/g^^j, in 



the CGHS mean-field theory (which, in particular, re- 
duces to A/g^^^j in the classical theory). This mass 
remains positive throughout the evaporation process of 
the mean-field approximation. 

• Although the horizon area goes to zero at the end of 
the evaporation process in the mean-field approxima- 
tion, A'/g^j'^j is not of Planck size at that time (i.e., at the 
point where the 'last ray' of Fig. 2 intersects 1^ ). For 
all black holes with large initial ADM mass, as the hori- 
zon area shrinks to zero A/g^^^j approaches a universal 
value sa 0.8647V in Planck units, with N = N/24. This 
end point Bondi mass is macroscopic since N is neces- 
sarily large in the semi-classical theory. 

• Dynamics during the evolution process also shows a 
universal behavior. For example, one can calculate 
-^Bondi ^ function of the horizon area ahor- There is 
a transient phase immediately after the horizon is first 
formed, though after that the plot of Afg^^j against 
ahor joins a universal curve all the way to zero area. 

• The flux of energy radiated across departs from the 
thermal flux when A/g^^^j and even ahor are macro- 
scopic. 

• Although the Ricci scalar of the mean-field metric g di- 
verges at the (weak) singularity, it is regular on the last 
ray and goes to zero as one approaches X^ along this 
ray. Thus, contrary to a wide spread belief, there is no 
'thunderbolt' curvature singularity in the semi-classical 
theory. 

We will see in section IV C that our results strongly sug- 
gest that the S matrix from to 1^ is hkely to be unitary. 
However, because of the universality of physical quantities at 
I^, it is very unlikely that information in the infalling matter 
at will be recovered in the outgoing state at X^. This is in 
sharp contrast with a wide-spread expectation; indeed, mech- 
anisms for information recovery have been suggested in the 
past (see e.g. [14]). This expectation illustrate the degree to 
which universality was unanticipated in much of the CGHS 
literature. 

The rest of this paper is organized as follows. In section 
II we summarize the mean-field theory of [7, 8], introduce 
the scaling behavior and explain the subtleties associated with 
Planck units in 2 dimensions. This framework provides the 
basis for the numerical results discussed in III. These results 
are more extensive than the brief summary presented above. 
Finally in section IV we discuss the ramifications of these re- 
sults and list interesting problems they suggest in geometric 
analysis. 



we will discuss in Sec. II B, it is actually the ratio A/adm/^ which classi- 
fies a solution, and more importantly there is qualitatively different behav- 
ior in the small (Planck size) vs. large (macroscopic) A/adm/^ Umits. In 
this terminology, the simulation of [12] corresponds to an initially Planck 
size black hole. 



n. THE MEAN-FIELD APPROXIMATION 

This section is divided into three parts. We recall (mainly 
from [2, 7, 8]) the PDEs that govern semi-classical gravity in 
the first part and their consequences in the second. In the third 
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we introduce a scaling symmetry which leads to a natural dis- 
tinction between macroscopic and Planck scale black holes. 



A. mean-field equations 

As mentioned in section I, in semi-classical gravity we have 
N quantum scalar fields with i = 1 . . . A^. In the mean- 
field approximation, we capture the idea that it is only the left 
moving modes that undergo gravitational collapse by choos- 
ing the initial state appropriately: we let this state be the 
vacuum state for the right moving modes /_ and a coher- 
ent state peaked at any given classical profile /° for each 

of the N left moving fields /|'''. This specification at 2^ 
defines a (Heisenberg) state j^"). Dynamical equations are 
obtained by taking expectation values of the quantum evo- 
lution equations for (Heisenberg) fields in this state I^I*) and 
ignoring quantum fluctuations of geometry but not of matter. 
Technically, this is accomplished by substituting polynomi- 
als P(8, <&) in the geometrical operators with polynomials 
P((9), ($)) := P($,9) of their expectation values. For the 
matter fields /^*\ on the other hand, one does not make this 
substitution; one keeps track of the quantum fluctuations of 
matter These lead to a conformal anomaly: While the trace 
of the stress-tensor of scalar fields vanishes in the classical 
theory due to conformal invariance, the expectation value of 
this trace now fails to vanish. Therefore equations of motion 
of the geometry acquire new source terms of quantum origin 
which modify its evolution. 

To summarize, then, in the mean-field approximation the 
dynamical objects are again just smooth fields /,;, 0, $ (rep- 
resenting expectation values of the corresponding quantum 
fields). While there are N matter fields, geometry is still en- 
coded in the two basic fields 9, $ which determine the space- 
time metric via g"'' = il?/"'' := Q^^^jf^. Dynamics 
of /('',9,$ are again governed by PDEs but, because of 
the trace anomaly, equations governing 9, $ acquire quan- 
tum corrections which encode the back reaction of quantum 
radiation on geometry. 

As in 4-dimensional general relativity, there are two sets 
of PDEs: Evolution equations and constraints which are pre- 
served in time. As in the classical theory, it is simplest to fix 
the gauge and write these equations using the advanced and re- 
tarded coordinates of the fiducial Minkowski metric. The 
evolution equations are given by 

□(r,)/''^'=0 ^ □(3)/«=0. (2.1) 

for matter fields and 

a+ 9_ $ + K^9 = G (f+_) = NGhd+ d- In $9^^2.2) 
$9+9-ln9 = -G(T+_) = -NGhd+ d-ln^Q-^ 

(2.3) 

for the geometrical fields where, as before, N ~ N/24:. The 
constraint equations tie the geometrical fields 9 , $ to the mat- 
ter fields fi. They are preserved in time. Therefore we can 



impose them just at T where they take the form: 

-dl $ + (9_ In 9 = G (f__) = (2.4) 

and 

-dl $ + 9+ $a+ In 9 = G (f++) = UNG {d+f°f 

(2.5) 

where = stands for 'equality atl^.' 

We will conclude this discussion of the field equations with 
a few remarks, and a description of our initial conditions. Be- 
cause /_ are all in the vacuum state, it follows immediately 
that, as in the classical theory, all the right moving fields van- 
ish; /I = also in the mean-field theory. Similarly, because 
Z^*'' are in a coherent state peaked at some classical profile 

/°, it follows that, for all i, /+\z+) = (on the en- 

tire fiducial Minkowski manifold M°). Thus, as far as matter 
fields are concerned, there is no difference between the clas- 
sical and mean-field theory. Similarly, as in the classical the- 
ory, we can integrate the constraint equations to obtain initial 
data on two null hypersurfaces. We will assume that /|°' van- 
ishes to the past of the line z+ = z+. Let denote the line 
z+ = z+ and the portion of the line z^=Zq"<C— l/^to 
the future of z+ = z+. We will specify initial data on these 
two surfaces. The solution to the constraint equations along 
these lines is not unique and, as in the classical theory we re- 
quire additional physical input to select one. We will again 
require that be in the dilaton vacuum to the past of /j^ and 
by continuity on I£ . Following the CGHS literature, we will 
take it to be $ = e'*^(^^-^ ). 5 Thus, the initial values of 
semi-classical 9, $ coincide with those of classical 9, $: 

Q^^^izt-z~) onallof /l and/p; (2.6) 

and 

$ = 9 on and, 

$ = 9 - 12NGjf^dz+ e^""^ /f ldl+ e"''^^^ ( 

on (2.7) 

(see (1.7)). The difference in the classical and semi-classical 
theories lies entirely in the evolution equations (2.2) and (2.3). 
In the classical theory, the right hand sides of these equations 
vanish whence one can easily integrate them. In the mean- 
field theory, this is not possible and one has to take recourse 
to numerical methods. Finally, while our analytical consider- 
ations hold for any regular profile /°, to begin with we will 
follow the CGHS Hterature in Sec. Ill A and Sec. IIIB and 
specify to represent a collapsing shell: 

12N (^] =MABuSiz+) (2.8) 



' Strictly, since <E> is an operator on tlie tensor product of N Fock spaces, one 
for eacli tlie expectation value is Ne""^-^^ \ But this difference 
can be compensated by shifting z~ . We have chosen to use the convention 
in the literature so as to make translation between our expressions and those 
in other papers easier. 
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so the shell is concentrated at z"*" = 0. In the literature this 
profile is often expressed, using in place of z+, as: 



127V 



(2.9) 



where {x+) = f'-°Hz+). In Sec. IIIC we will discuss 
results from a class of smooth matter profiles. 



B. Singularity, horizons and ttie Bondi mass 

The classical solution (1.7) has a singularity £s where <& 
vanishes. As remarked in section I, in the mean-field theory, 
a singularity persists but it is shifted to $ = 2NGh [2]. The 
metric g"'' = Q^^^rj"'^ is invertible and continuous there but 
not C^. Thus the singularity is weakened relative to the clas- 
sical theory. Furthermore, its spatial extension is diminished. 
As indicated in Fig. 2, the singularity now originates at a finite 
point on the collapsing shell (i.e. does not extend to X^) and 
it ends in the space-time interior (i.e., does not extend to ). 

What is the situation with horizons? Recall from sec- 
tion I that, in the spherically symmetric reduction from 4- 
dimensions, = e^^'^/k^ := and each round 2-sphere 
in 4-dimensional space-time projects down to a single point on 
the 2-manifold M. Thus, in the CGHS model we can think of 
$ as defining the 'area' associated with any point. (It is di- 
mensionless because in D space-time dimensions the area of 
spatial spheres has dimensions [L]^^"^.) Therefore it is natu- 
ral to define a notion of trapped points: A point in the CGHS 
space-time {M, g) is said to he future trapped if 9+$ and 9-$ 
are both negative there and future marginally trapped if 
vanishes and c)_$ is negative there [2, 15]. In the classical 
solution resulting from the collapse of a shell (2.8), all the 
marginally trapped points lie on the event horizon and their 
area is a constant; we only encounter an isolated horizon [16] 
(see Fig.l). The mean-field theory is much richer because it 
incorporates the back reaction of quantum radiation. In the 
case of a shell collapse, the field equations now imply that a 
marginally trapped point first forms at a point on the shell and 
has area 



ainitial (* - 27VG?l) I initial 



-NGh + NGh 1 



^"ADM 



(2.10) 



As time evolves, this area shrinks because of quantum radi- 
ation [2]. The world-line of these marginally trapped points 
forms a generalized dynamical horizon (GDH), 'generalized' 
because the world-line is time-like rather than space-like [16]. 
(In 4 dimensions these are called marginally trapped tubes 
[17].) The area finally shrinks to zero. This is the point at 
which the GDH meets the end-point of the (weak) singularity 
[2, 12, 18] (see Fig. 2). It is remarkable that all these inter- 
esting dynamics occur simply because, unlike in the classical 
theory, the right sides of the dynamical equations (2.2), (2.3) 
are non-zero, given by the trace-anomaly. 



We will see in section III that while the solution is indeed 
asymptotically flat at 1^ , in contrast to the classical solution, 
Zj^ is no longer complete. More precisely, the space-time 
(Af , g) now has a future boundary at the last ray — the null 
line to from the point at which the singularity ends — and 
the affine parameter along with respect to gab has a finite 
value at the point where the last ray meets X^. Therefore, 
in the semi-classical theory, we cannot even ask if this space- 
time admits an event horizon. While the notion of an event 
horizon is global and teleological, the notion of trapped sur- 
faces and GDHs is quasi-local. As we have just argued, these 
continue to be meaningful in the semi-classical theory. What 
forms and evaporates is the GDH. 

Next, let us discuss the structure at null infinity [7, 8]. As in 
the classical theory, we assume that the semi-classical space- 
time is asymptotically flat at X^ in the sense that, as one takes 
the limit z+ — >^ cx) along the lines — const, the fields $, Q 
have the following behavior: 



$ = A{z~) e«^^ + B{z-) + Oie-"'^) 
e = A(z-) e«^* + B(z-) + 0(e-''^^) , 



(2.11) 



where A, i3,A, B are some smooth functions of z~. Note 
that the leading order behavior in (2.11) is the same as that 
in the classical solution. The only difference is that i?, B are 
not required to be constant along because, in contrast to 
its classical counterpart, the semi-classical space-time is non- 
stationary near null infinity due to quantum radiation. There- 
fore, as in the classical theory, X^ can be obtained by taking 
the limit z+ — )• oo along the lines z~ = const. The asymp- 
totic conditions (2.1 1) on 8, $ imply that curvature — i.e., the 
Ricci scalar of gab — goes to zero at X^ . We will see in section 
III that these conditions are indeed satisfied in semi-classical 
space-times that result from collapse of matter from X^ . 

Given this asymptotic fall-off, the field equations determine 
A and B in terms of A and B. The metric gab admits an 
asymptotic time translation which is unique up to a con- 
stant rescaling and the rescaling freedom can be eliminated by 
requiring that it be (asymptotically) unit. The function A{z~) 
determines the affine parameter of t° via: 



-Ky 



A{z- 



(2.12) 



Thus y~ can be regarded as the unique asymptotic time pa- 
rameter with respect to gab (up to an additive constant). Near 
the mean-field metric g can be expanded as: 



dS-^ = - (l + Be'^^y'-y^ + Oie-^'^y )j dy+ dy- 

(2.13) 

where j/+ = z+. 

Finally, equations of the mean-field theory imply [7, 8] that 
there is a balance law at X^: 



' dz 2 dz 

NhG rdV ]2 
2 ^dz-^^dz-' J 



(2.14) 
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In [7], this balance law was used to introduce a new notion 
of Bondi mass and flux. The left side of (2.14) led to the 
definition of the Bondi mass: 



ATV 
Bondi 



ay dz ^ az ■ 



while the right side provided the Bondi flux: 



7ATV 



so that we have: 



NhG rd'^y- dy- . 
2 ''dl^' 



ATV 
Bondi 



dy 



ATV 



(2.15) 



(2.16) 



(2.17) 



By construction, as in 4 dimensions, the flux is manifestly 
positive so that -A/g^-,^; decreases in time. Furthermore, it 
vanishes on an open region if and only if y^ = Ciz^ + C2 
for some constants Ci, C2, i.e. if and only if the asymptotic 
time translations defined by the physical, mean field metric g 
and by the fiducial metric r/ agree at 1^ , or, equivalently, ;/ 
and only if the asymptotic time translations of g on and 
agree. Finally, note that g"'' = /± = 0, $ = 9 = 
exp K,{z^ — z~), is a solution to the full mean-field equations. 
As one would expect, both A/g^^^j and F^^^ vanish for this 
solution. 

The balance law is just a statement of conservation of en- 
ergy. As one would expect, h appears as an overall multi- 
plicative constant in (2.16); in the classical theory, there is no 
flux of energy at . If we set h = Q, M^J^^:^ reduces to the 
standard Bondi mass formula in the classical theory (see e.g., 
[15]). Previous literature [1, 2, 14, 15, 18-20] on the CGHS 
model used this classical expression also in the semi-classical 
theory. Thus, in the notation we have introduced here, the 
traditional definitions of mass and flux are given by 



Trad 
Bondi 



dB 
dy- 



kB, 



and 



^Trad _ ^ATV 



NhG 



d ^ d^y ^dy ^^2- 



dy dz 



'dz- 



(2.18) 



(2.19) 



As noted in section I, numerical simulations have shown that 
^^Bondi become negative and large even when the horizon 
area is large, while M^J^^-^ remains positive throughout the 
evaporation process. 



C. Scaling and the Planck regime 

Finally, we note a scaling property of the mean-field theory, 
which Ori recently and independently also uncovered [11] and 
which is also observed in other quantum gravitational sys- 
tems [21]. We were led to it while attempting to interpret 
numerical results which at first seemed very puzzling; it is 
thus a concrete example of the how useful the interplay be- 
tween numerical and analytical studies can be. Let us fix z^ 



and regard all fields as functions of z^. Consider any solu- 
tion (0, $, N, f^^^) to our field equations, satisfying bound- 
ary conditions (2.6) and (2.7). Then, given a positive number 
A, (e,$,iV,/f )givenby'' 

InA. 



A$(z^ 



+ ), N = \N 



K 

In A, 



-), 



is also a solution satisfying our boundary conditions, where, 
as before, we have assumed that all scalar fields have an iden- 
tical profile Note that /° is completely general; we have 
not restricted ourselves, e.g., to shells. Under this transforma- 
tion, we have 



ab 



ah 



■In A 



-^^ADM — > AA/aDM 

'^-''^ Bondi 
A^^ATV 



ywrATV 
-'^^ Bondi 
^ATV 

acDH 



(2.20) 



where acDH denotes the area of the generalized dynamical 
horizon. This symmetry implies that, as far as space-time ge- 
ometry and energetics are concerned, only the ratios M /N 
matter, not separate values of M and N themselves (where 
A/ can either be the ADM or the Bondi mass). Thus, for ex- 
ample, whether for the evaporation process a black hole is 
'macroscopic' or 'Planck size' depends on the ratios AI/N 
and a.Guti/N rather than on the values of M or acDH them- 
selves. 
We will set 



M* = Maum/N 
il^Bondi = MilnJN, and 
m = ^^Bondillast ray 



(2.21) 



(We use N = N/24: in these definitions because the dynami- 
cal equations feature N rather than N.) We will need to com- 
pare these quantities with the Planck mass. Now, in 2 dimen- 
sions, G, h and c do not suffice to determine Planck mass, 
Planck length and Planck time uniquely because Gh is dimen- 
sionless. But in 4 dimensions we have unambiguous defini- 
tions of these quantities and, conceptually, we can regard the 
2 dimensional theory as obtained by its spherical reduction. In 
4-dimensions, (using the c=l units used here) the Planck mass 
is given by A/pj = /1/G4 and the Planck time by Tpj = G4 h. 
From Eqs (1.2) and (1.3) it follows that G4 is related to the 2- 
dimensional Newton's constant G via G = GiU^. Therefore 
we are led to set 



AT, 



PI 



~G' 



and Tpj = — — 



Gh 

2" 



K 



(2.22) 



The shift in 
gKj^j,+ — z") on 
See footnote 5. 



is needed because we chose to use the initial value © = 
as in the literature rather than = Ne"^^^ ~^ \ 
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When can we say that a black hole is macroscopic? One's 
first instinct would be to say that the ADM mass should be 
much larger than Mp\ in (2.22). But this is not adequate for 
the evaporation process because the process depends also on 
the number of fields N. In the external field approximation 
where one ignores the back reaction, we know that at late 
times the black hole radiates as a black body at a fixed tem- 
perature Tnaw = nh. ^ The Hawking energy flux at is 
given by = NK'^h/2. Therefore the evaporation pro- 

cess will last much longer than 1 Planck time if and only if 
{Madm/F^'^'^) > Tpi, or, equivalent^ 

M* > Gh Mpi. (2.23) 

(Recall that Gh is the Planck number) So, a necessary con- 
dition for a black hole to be macroscopic is that M* should 
satisfy this inequality. In section III we will see that, in the 
mean-field theory, quantum evaporation reveals universality 
ab-eady if M* >AGhMpi. 



III. ANTICIPATED AND UNFORESEEN BEHAVIOR 

All physical predictions of the mean-field theory arise from 
the set of 5 equations (2.1) - (2.5). The only difference from 
the classical theory lies in the fact that, because of the trace 
anomaly, right hand sides of the dynamical equations (2.2) 
and (2.3) are no longer zero. But this difference has very sig- 
nificant ramifications. In particular, it is no longer possible to 
obtain explicit analytical solutions; one has to take recourse to 
numerics.** 

Also, we now have interesting time-dependent phenomena 
such as formation and evaporation of dynamical horizons. 
Given these differences, a number of global questions natu- 
rally arise. Does the space-time continue to be asymptotically 
flat at in the mean-field theory? If so, is complete as 
in the classical theory? Is the Bondi mass positive? Does it 
go to the Planck scale as the horizon area goes to zero? Is 
the flux of energy of the quantum radiation constant, given by 
hn^N /2 at late times, as in the external field approximation 
a la Hawking [25, 26]? If not, the quantum radiation would 



Note that this relation is the same as that in 4 dimensions because the clas- 
sical CGHS black hole is stationaiy to the future of the collapsing matter 
with surface gravity k. However, there is also a key difference: now k 
is just a constant, independent of the mass of the black hole. Therefore, 
unlike in 4-dimensions, the temperature of the CGHS black hole is a uni- 
versal constant in the external field approximation. Therefore, when back 
reaction is included, one does not expect a CGHS black hole to get hotter 
as it shrinks. 

There are variants of the CGHS model that are explicitly soluble, for exam- 
ple the RST (Russo-Susskind-Thorlacius) model [22], and the Bilal-Callan 
model [23]. However, results obtained in these models are not likely to be 
generic even in 2 dimensions because of their extra symmetries [12, 20]. 
More importantly, it was pointed out in [2, 12, 24] that the RST model 
is inconsistent even in the large A'^ limit, and the Bilal-Callan model has 
a Hamiltonian that is unbounded from below. Thus though they exhibit 
many features of general 2D semi-classical black hole evaporation, they 
are physically less interesting than the CGHS model. 



not be compatible with thermal flux, violating a key assump- 
tion that has been made over the years. 

Recall from section II B that since our profile functions 
vanish to the past of a null ray , the solution in the past is 

given by = 0, gab = Vab, 6 = $ = e^<^'^-''\ To 
obtain it to the future of , we use the initial data on and 
given by (2.6) and (2.7). The dynamical equations (2.1), 
(2.2) are hyperbolic and therefore well adapted to numerical 
evolution. They have been solved numerically before (see, 
in particular, [12, 18-20]), and much information has been 
learnt about the CGHS model, for example the dynamics of 
the GDH and that it evaporates to zero area, terminating in 
a naked singularity. However, to our knowledge in all these 
simulations, the choice of parameters A/adm and N was such 
that M* = Madm/N was less that 2.5. As we will see, in a 
precise sense, these black holes are Planck scale already when 
they are formed. 

Reliable simulations of macroscopic black holes with 
?i 5 — 10 turn out to be much harder to perform, and 
several additional steps beyond a straight-forward discretiza- 
tion are necessary to study this regime [9]. First, one needs 
to formulate the problem in terms of 'regularized' variables 
which do not diverge at infinity. Second, one needs to intro- 
duce coordinates which (a) bring the infinite portion of 
of interest to a compact interval, and (b) enlarge the region 
near the last ray by a factor of roughly e^^ relative to a uni- 
form discretization of the (compactified) time-translation co- 
ordinate z". (In other words, a region in the vacuum solution 
nearly of physical size equal to ^ e~^' expands to a region 
of physical size 0(1) on I^, where all the interesting dynam- 
ics occurs.) Third, to achieve this latter part of the coordinate 
transformation one needs to know the location of the last ray 
extremely well, requiring high accuracy numerical methods. 
This is achieved by using Richardson extrapolation with in- 
termittent error removal, beginning with a unigrid method that 
is second order accurate. With four successively finer meshes 
the overall rate of convergence of the scheme is 0{h'^), and 
this was sufficient to reduce the truncation error to the order 
of machine round-off 10^^^). These high precision nu- 
merical calculations showed that, while some of the long held 
assumptions on the nature of quantum evaporation are borne 
out, several are not. We summarize these results from a physi- 
cal perspective in the next three sub-sections. Complementary 
discussion of numerical issues appears in [9]. 

Numerical calculations were performed for a range of 
rescaled ADIVI masses M* from 2^^° to 16 and TV varying 
from 12 to 24000. Since in this paper we are primarily inter- 
ested in black holes which are initially macroscopic, we will 
focus on A/* > 1 and, since our simulations exhibited the 
expected scaling behavior, all our plots, except Fig. 5, are for 
= 24 (i.e., N = 1). Finally, in these simulations we set 
h ^ G = K = 1. The first two sub-sections summarize re- 
sults from a shell collapse and the third sub-section reports on 
results obtained from more general infalling profiles. 
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8 12 16 20 

Figure 3. The Ricci scalar R for Af* = 8. Left: 2D contour plot 
of R^^^ showing the increase in R as the singularity (dark vertical 
region near the middle) is approached and the asymptotically flat re- 
gion (7? — 0) near — >■ oo). Right: R^^^ as a function of 
z+ on the lines z~ - aT^^ = -10"'^, -10"**, -10"* (marked on 
the left panel as horizontal lines), showing a double peak, indicating 
the divergent behavior of d+d-^ at the singularity. The fact that 
the peak is narrow rules out a strong thunderbolt singularity. Note 
that the dark color at the region of the singularity is due to the high 
density of contour lines, and not directly due to negative values of 
R. While naive numerical calculation of R very close to does 
not yield reliable results due to catastrophic cancelation, it is already 
very small in the high z"*" values shown here, and the trend towards 
is clear. 



A. Shell Collapse: Anticipated Behavior 

Asymptotic flatness at I^: First, 0,<i> do indeed satisfy 
the asymptotic conditions (2.11). This was also noted in 
the recent approximate solution to the CGHS equations by 
Ori [11]. The simulations provide values of the functions 
A{z~)^B{z~) and y^{z^). As a consistency check on the 
simulation, we verified the balance law (2.14) by calculating 
separately the right and left sides of this equation as close to 
the last ray as the numerical solution gave reliable (conver- 
gent) results. We also computed the scalar curvature R of the 
mean-field metric g, and it does go to zero at — see Fig. 3 
for an example. 

Finiteness of at the last my: In the classical solution, 
the affine parameters y~ along and z~ along defined 
by g are related by 



GM 



(3.1) 



Hence y 



oo at Kz 



\n{GM / k). This is the point 



at which the singularity and the event horizon meet (see 
Fig 1). Thus, in the classical solution is complete but, 
in a precise sense, smaller than . For a test quantum field 
/_ on the classical solution, one then has to trace over modes 



N 
T3 



T3 



O 



< 



N 

1^ 



-10 



^10 




^10 



Figure 4. Left: Plot of logjo (dy /Az ) vs logjo A for M* = 8, 
iV = 1, where A — (z^^^^g — 2^). Right: Slope of the curve on 
the left. If locally the function on the left behaves as ^ (kA)~^, the 
curve on the right shows —p. In the distant past (rightmost region in 
both plots), y~ tends to z~ . The intermediate region is similar to that 
in the classical solution where {dy~ /dz~) ~ (kA)~^. As the last 
ray is further approached (leftmost region), we see that {dy~ /dz~) 
increases much slower than («:A)~^, leading to a finite value for y~ 
at the last ray. 



residing on the part of Ij^ which is missing from I^. This 
fact is directly responsible for pure states on X£ to evolve to 
mixed states on 1^ , i.e., for the non-unitarity of the S'-matrix 
[7, 8] of the test field. What is the situation in the mean-field 
theory? Our analysis shows that, as generally expected, the 
affine parameter w.r.t. the mean field metric g takes a finite 
value at the last ray on ; a necessary condition for unitarity 
of the S -matrix is met. 

To establish this result, we apply the following strategy. Let 
us return to the classical solution g and set 

'*^sT„g,ci = - HGM/ n) and = zT^^^^j - z" . 

(3.2) 

(The subscript 'sing,cr just highlights the fact that this is the 
point at which the classical singularity meets .) Then we 
have 



= - - In (1 



(3.3) 



When Aci tends to zero, is dominated by the leading or- 
der term — (1/k) 1ii(— kAci) which diverges at Ad = 0. This 
logarithmic divergence is coded in the power —1 in the ex- 
pression of the derivative (dy^/dz~): 



^ = (AjAel) 

dz 



finite terms 



(3.4) 



If we had (k Ad)"^ on the right side rather than (k Ad)^^, 
then y~ would have been finite at the future end of of g 
forp < 1 (as then y' = (k Ad)^~^/(1 — p)+ finite terms). 
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In the mean-field theory, the last ray starts at the end point 
of the singularity and meets of g at its future end point. 
We will denote it by the line z~ = •2:^ng- Following the above 
discussion, to show that the affine parameter w.r.t. g is fi- 
nite at z~ = 2g[jjg we focus on the behavior of (d?/~/dz~) 
near this future end point of . More precisely, we analyze 
the functional behavior of (dj/~/dz~) and determine a lo- 
cal p extracted from the logarithmic derivative of {dy^ / dz^) 
with respect to A = z^T^^ — . Results in Fig. 4 show that 
{dy~ /dz~) grows much slower near the last ray in the mean- 
field theory than it does in the classical theory. In fact, over the 
entire range of the local estimate of p is strictly less that 
1, and asymptotes to approaching the last ray. This implies 
that is finite at the last ray in the mean-field theory. 

Note that the above analysis is only valid if we have de- 
termined the location of the singularity with sufficient accu- 
racy such that the numerical uncertainty in its location is much 
smaller than the range in A where we extract the asymptotic 
behavior of the function. From convergence studies we es- 
timate our precision in determining zT^^^ is at the order of 
10"^'^, and hence all the values in Fig. 4 are sufficiently far 
from the last ray to provide a reliable measure of the power p. 



B. Shell collapse: Unforeseen Behavior 

The numerical calculations also revealed a number of sur- 
prises which we now discuss. 

Botidi mass for large N: Scaling properties discussed in 
section 11 imply that if the Bondi mass at the last ray is finite, 
it will be macroscopic for a sufficiently large N. This expec- 
tation is borne out (in particular the Bondi mass is finite) in 
all our simulations with large A/adm and large N. Fig. 5 
summarizes the result of a simulation where N — 720 and 
Madm = 360 (so iV = 30 and M* 12). The Bondi 
mass, MgQ^j^j, that has been commonly used in the litera- 
ture [1, 2, 14, 15, 18-20] becomes negative even far from the 
last ray when the horizon area is still macroscopic, and has 
a macroscopic negative value at the last ray.** On the other 
hand, the more recent M^J^^^ [7, 8] remains strictly positive. 
As one would expect from the scaling relations, because N is 
large, Mg^^j is also macroscopic at the last ray. 

Universality of the end state: Fig 6 shows a plot of m*, 
the value of {M^^^^JN) at the last ray, against M* = 
{Madm/ N), for several values of the initial M* > 1. The 
curve that fits the data, shown in the figure, is 



a {I 



') 



(3.5) 



with specific values for the constants a, /?, 7 

aw 0.864, /? w 1.42, 7^ 1.15. 



' After this work was completed, Javad Taghizadeh Firouzjaee pointed out 
to us that the fact that the traditional Bondi mass can become negative was 
already noticed in [19]. Again though, in our terminology the numeri- 
cal simulation in that work corresponds to a microscopic black hole with 
M* = IMpi. 




Sing 



GDH 



Figure 5. The ATV Bondi mass Afg^^j (solid lines) and the 
traditional Bondi mass M^^^^ (dashed lines) are plotted against 
z~ — z^ng (left) and the horizon area (right). This simulation corre- 
sponds to Madm = 360, = 720 (so M* = 12). For high values 
of A^, both formulas give a large non-zero Bondi mass at the last ray. 
-^^Bondi becomes negative when the area is still macroscopic. On the 
other hand Afg^^; remains strictly positive all the way to the last 
ray, where the generalized dynamical horizon (GDH) shrinks to zero 
area. 
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Figure 6. The value of m* (i.e. M-^J2di/^ ths l^st ray) is plotted 
against M* (which equals Madm/N) for M* > 1. For Macro- 
scopic Af* (actually, already for M* > 4!) m* has a universal 
value, approximately 0.864. 



It is visually clear from the plot that there is a qualitative dif- 
ference between M* > 4 and M* < 4. Physically this can be 
understood in terms of ainitiah the area of the first marginally 
trapped surface: Eq (2.10) implies that a*jjj^jg^[ = ajnitiai/^ 
can be greater than a Planck unit only if M* is larger than 
3. It is therefore not surprising that Af* = 4 should serve as 
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Figure 7. Left: The affine parameter (defined in Eq. (2.12) 
of the physical metric g is plotted against the rescaled area a* := 
(acDH / N) of the generalized dynamical horizon (given by ($/iV — 
2)) at the horizon for values of M* from 4 to 14. Even though the 
curves are very similar in shape, they do not coincide. Right: Once 
the shifting freedom in j/~ is utilized, we see that a properly shifted 
version is universal with respect to a* for all macroscopic M* 
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Figures. F* = (F^^^/iV) is plotted against the horizon area a* := 
[slgok/N) (left) and (right) for values of M* from 4 to 14. For 
all M* values, F* starts at the value of at the distant past (ki/^ ^ 
— 1), and then joins a universal curve of F* . Note that once the GDH 
is formed, (the rightmost beginning of each curve on the left plot) 
F* is already slightly larger in magnitude than the Hawking/thermal 
value 0.5 and it increases steadily as one approaches the last ray (i.e. 
as acDH and j/^^i approach 0. 



the boundary between macro and Planck regimes. Indeed, as 
Fig 6 shows, if M* > 4, the value of the end point Bondi 
mass is universal, m* w 0.864. For M* < 4 on the other 
hand, the value of m* depends sensitively on A/*. This could 
have been anticipated because if A/* < 3, what evaporates is 
a GDH which starts out with one Planck unit or less of area 
a"*. Thus, in the mean-field approximation it is natural to re- 
gard CGHS black holes with M* > A as macroscopic and 
those with M* < 4 microscopic. Finally, for macroscopic 
black holes, the end-value of the traditional Bondi-mass is 
also universal: A/T^^^. < ^^^^ {MZfjN) ^ -2.0 
as ahor 0. 

As noted in the beginning of section III, there have been 
a number of previous numerical studies of the CGHS model 
[12, 18-20]. They clarified several important dynamical is- 
sues. However they could not unravel universality because 
they all focused on cases where the black hole is microscopic 
already at its inception: M* < 2.5 in [18], M* = 1 in [12] 
and [19] and M* =- 0.72 in [20]. This limitation was not no- 
ticed because the scaling symmetry and its significance was 
not appreciated. 

Dynamical universality of y^: The horizon area aooH 
(more precisely, its negative) provides an invariantly defined 
time coordinate to test dynamical universality of other phys- 
ical quantities. Let us begin with , the affine parame- 
ter along with respect to the physical metric g defined 
in (2.12). Fig. 7, left, shows the plot of y^ against a* := 
(slg-dh/N) for various values of M*. These plots show that 



the time dependence of y~ for various values of M* is very 
similar but not identical. Recall, however, that there is some 
freedom in the definition of the affine parameter In particu- 
lar, in each space-time we can shift it by a constant, and the 
particular value of the constant can vary from one space-time 
to the next (e.g. depend on the ADM mass). This shift does 
not affect any of our considerations, including the balance law 
(2.14). 

Let us define y^^^ by shifting each y^ so that each solution 
reaches the same small non-zero value of the horizon area, 
a"^ = e, at the same y^. It turns out that this shift has the 
remarkable feature that, for initially macroscopic black holes, 
all shifted curves now coincide for all values of a*. Thus, we 
have a universal, monotonic function of a* plotted in Fig. 7, 
right. Hence y^^^ also serves as an invariant time coordinate. 
In fact it has an advantage over acoH: whereas a* is defined 
only after the first marginally trapped surface is formed (see 
Fig 2), y~^ is well defined throughout the mean-field space- 
time (Af, g). 

Dynamical Universality of F^^^ and A/g^^^j: We can re- 
peat the procedure used above for y" to investigate if dynam- 
ics of other physical quantities such as the Bondi flux F* := 
(F^TV/^) and the Bondi mass M^^^^^, := (M^JZi/N) 
are also universal. Note, however, that unlike y , there is 
no 'shift' (or indeed any other) freedom in the definitions 
of F^^^ and M^J^^. So, if there is universality, it should 
emerge directly, without any adjustments, in the plots of F* 
and ^-fBondi against a* = (acDn/^) or y^^^. 

Let us begin with the Bondi flux. Recall, first, that in the ex- 
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Figure 9. M^^^^-^ = (M-^Jj^-JN) is plotted against ttie horizon 
area a* ~ (acDnZ-^V) (left) and j/^ (right) for values of M* from 
4 to 14. For all these macroscopic M*, A/bq^^j; starts at the value of 
Madm in the distant past (kj/^^ <C and then joins a universal 
curve of M^^^^^^. When the dynamical horizon first forms A/gondi 
is quite close to its initial value of A/*, (This is difficult to see in the 
left plot where all the curves crowd.) This means that almost all of 
the evaporation occurs after the formation of the dynamical horizon. 

ternal field approximation [2, 26], the energy flux is very small 
in the distant past, rises steeply at Ky~ w — \n{GMABM/ 
and then quickly asymptotes to the Hawking value F^aw _ 
{NhK^/2). This constant flux is characteristic of thermal ra- 
diation at temperature kH in two space-time dimensions. In 
our simulations (with N = 24, or) TV = 1 and h = k = 1, it 
corresponds to i^Haw ^ q g 

In the mean-field theory, numerical simulations show that, 
for all initially macroscopic black holes, the energy flux 
F* := {F^"^^ / N) is also negligibly small in the distant past 
and then rises steeply. But this rise is now associated with a 
clearly identifiable dynamical process: formation of the first 
marginally trapped surface. As we noted in section 11 B, for 
a shell collapse, analytical calculations show that the area of 
this first surface is given by (2.10). Assuming that we have 
a macroscopic initial mass, M* ^ V Gh A/pi =: A/pL, Eq 
(2.10) simplifies: 

aL.ai«GM£-l + ^ + ...] (3.6) 

This relation is borne out in simulations. Assuming that the 
black hole is very large at this stage, heuristically, one can 
equate the area of this new born GDH with the Bondi mass 
at the retarded instant of time, say y~ = y^, at which is it 
born. This implies that, per scalar field, only ~ 1 Planck unit 
of A/gjjj^jj has been radiated over the long period of time from 
y" = — oo till y^ ~ Vo ■ ^^^t once the GDH appears, the flux 
rises steeply to a value close to but higher than 0.5. Then, it 
joins a universal curve all the way to the last ray where the area 



a* shrinks to zero. Thus, after a brief transient phase around 
the time the GDH is first formed, the time-dependence of the 
Bondi flux is universal. Fig. 8, left shows this universal time 
dependence with a* as time and Fig. 8, right shows it with y~^ 
as time. 

In virtue of the balance law (2.14) one would expect this 
universality to imply a universal time dependence also for the 
Bondi mass A/g^j^^^j. This is indeed the case. At spatial in- 
finity ij^, we have Afg^j^^j = A/*. There is a transient phase 
around the birth of the GDH in which the Bondi mass de- 
creases steeply. Quickly after that, the time dependence of 
^^Bondi follows a Universal trajectory until the last ray. Fig. 9, 
left shows this universality with a* as time while Fig. 9, right 
shows it with y~^ as time. 

To summarize, using either a* or y^j^ as an invariant time 
coordinate, we can track the dynamics of F* and A/g^j^^j. In 
each of the four cases, there is a universal curve describing 
these dynamics. For definiteness let us use a* as time and 
focus on A/gQjjj; (the situation is the same in the other three 
cases). Since both quantities are positive, let us consider the 
time-mass quadrant they span. Fix a very large initial black 
hole with A/* = A/* and denote by Cq the curve in the quad- 
rant that describes its time evolution. Then, given any other 
black hole with M* < M*, the curve c describing the dynam- 
ical evolution of its Afgonj; starts out at a smaller value of time 
(i.e. a*) marking the birth of the GDH of that space-time and, 
after a brief transient phase, joints the curve Cq all the way 
until its horizon shrinks to zero. Here we have focused on 
the ATV flux and mass because their properties make them 
physically more relevant. But this universality holds also for 
the flux and mass expressions, i^Ti'^^^ -^Bondi '^hat have been 
traditionally used in the literature. 

Curvature at the last ray: There has been considerable dis- 
cussion on the nature of geometry at the last ray. Since this 
ray starts out at the singularity, a natural question is whether 
a curvature singularity propagates out all along the last ray 
to . This would be a 'thunderbolt' representing a singular 
Cauchy horizon [20]. If it were to occur, the evolution across 
the last ray would not just be ambiguous; it would be impossi- 
ble. However, a priori it is not clear that a thunderbolt would 
in fact occur. For, the 'strength' of the singularity goes to zero 
at its right end point where the last ray originates. 

Using numerical simulations. Hawking and Stewart [20] ar- 
gued that a thunderbolt does occur in the semi-classical the- 
ory. But they went on to suggest that it could be softened 
by full quantum gravity, i.e., that full quantum gravity effects 
would tame it to produce possibly a very intense but finite 
burst of high energy particles in the full theory. 

Our calculation of the Ricci scalar very close to the last 
ray shows that, except for a small region near the singularity, 
the scalar curvature at the last ray is not large (Fig. 3). Thus, 
our more exhaustive and high precision calculations rule out 
a thunderbolt singularity in the original sense of the term. 
This overall conclusion agrees with the later results in [18]. 
(Both these calculations were done only for initially micro- 
scopic black holes while results hold also for macroscopic 
ones.) However, our calculations show that the flux F^^^ 
does increase very steeply near the last ray (see Fig. 8). Nu- 
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merically, we could not conclude whether the flux remains 
finite at the last ray or diverges. However, the integrated flux 
which determines the change in A/b^^i is indeed finite and 
in fact not very significant in the region very near the last ray. 
For macroscopic M* values, the total radiated energy after the 
point when F* reaches the value 1 is ~ 1 Planck mass, (see 
Figs. 8, 9). Thus, if we were to associate the thunderbolt idea 
to the steep increase of flux at the last ray, this would have 
to be in quite a weak sense; in particular, there is no singular 
Cauchy horizon. 

Nature of the Bondi flux: Recall that in the external 
field approximation, the energy flux starts out very low, 
rapidly increases and approaches the constant thermal value 
{F^^™ /N) = (= 0.5 in our simulations) from below 

[2, 26]. In the mean-field theory, the flux F^'^^ also starts 
out very small and suddenly increases when the GDH is first 
formed. However, it overshoots the thermal value and ceases 
to be constant much before the black hole shrinks to Planck 
size (Fig. 8). During subsequent evolution, F^^^ monotoni- 
cally increases in magnitude and is about 70% greater than the 
constant thermal value F^^'" when Memidi ^ 2NMpi: the 
standard assumption that the flux is thermal till the black hole 
shrinks to Planck size is not borne out in the mean field theory. 
(One's 4-dimensional intuition may lead one to think that the 
increase in the flux merely reflects that the black hole gets hot- 
ter as it evaporates; but this is not so because the temperature 
of a CGHS black hole is an absolute constant, Tnaw = kH). In 
the interval between the formation of the GDH and the time 
when Mg^^j approaches TVMpi, the numerical flux is well 
approximated by 
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Figure 10. The value of m* (i.e. M-^J^^^^/N at the last ray) plot- 
ted against M* (which equals Maom/ N) for M* > 1. In addition 
to points corresponding to shell collapse (w — 0) the plot now in- 
cludes points with more general profiles with w — 0.25, 0.5, 1. The 
universal value m* ~ 0.864 persists for M* > 4. 



test this conjecture, we evolved a 2-parameter family of ini- 
tial data, parameterized by a characteristic initial mass M and 
width w. Now, it is clear from the form (2.6), (2.7) of initial 
data that what matters is not the profile itself but rather 

the integral of (9+/^°^)^. We will specify it using two param- 
eters, M and w: 
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Thus, in this interval the flux is close to the constant thermal 
value only while the area a of the GDH is much greater than 
N Planck units.'" 



C. Universality beyond the shell collapse. 

So far, we have focused our attention on a delta distribution 
shell collapse (2.8). As we will discuss more in the following 
section IV, we expect the results to be robust for a large class 
of infalling profiles, so long as the GDH forms promptly. To 



'"The leading order correction + {N Mpi/ M^^^^^) to the Hawking flux 
was obtained by Ori by analytical approximation methods and served as 
the point of departure for obtaining the fit (3.7). Note also that if the fluxes 
differ over a significant time interval, it follows that the quantum radiation 
is not thermal. But the converse is not true as there are pure states in the 
outgoing Hilbert space for which the energy flux at is extremely well 
approximated by the constant thermal value. For quantum states, what 
matters is the comparison between the function y~i-^{z^) and its classical 
counterpart given by (3.1) [7, 8]. And these two functions are 

very different. Finally, nonthermal fluxes were also observed in a quantum 
model of four-dimensional spherical shell collapse [27] 
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(3.8) 

This choice is motivated by the following considerations. 
First, as in the shell collapse, there is a neighborhood of in 
which the solution represents the vacuum of the theory. Sec- 
ond, the power 4 on the right side is chosen to ensure high 



0). Thus, f^' is and 



differentiability at x"*" = 1 (i.e. z 

furthermore decays faster than e~*-^"^^ for any (7 as z+ — > oo. 
Third, the parameter w provides a measure of the width of the 
matter profile in coordinates, which is roughly the width 
in z+ for w < 1. Finally, note that we recover the shell profile 
in the limit w — and expect that the physical requirement 
of a 'prompt collapse' will be met for sufficiently small w. In 
the case of a shell profile (2.8), the parameter M represents 
the ADM mass. A simple calculation shows that for a general 
profile in family (3.8), Maum is given by a function of the 
two parameters: Madm = M{1 + 1.39 w). Thus, within this 
family, the issue of universality of a physical quantity boils 
down to the question of whether it depends only on the spe- 
cific combination + 1.39 w) of the two parameters. 

Numerical evolutions were carried out for M* sa 
6, 9, 11, 13 and w ~ 0.25, 0.5, 1. We find that universality is 
indeed retained for all these cases. Specifically, we repeated 
the following analysis of section IIIB for various values of 
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Figure 11. _F* (left) and Meondi (right) plotted against y^^^, for var- 
ious incoming matter profiles (w and Madm values), including sev- 
eral shell (w — 0) cases. The time when the dynamical horizon first 
forms is marked on each flux curve (which is later for larger w). All 
the curves with the same A/adm (6 in this example) are on top of 
each other and cannot be distinguished by the eye, showing that they 
have the same universal behavior throughout the evolution, including 
the early times. More generally all the asymptotic physical quantities 
depend only on the combination A/adm of the profile parameters A4 
and w as long as nw is small compared to the initial area of the GDH. 



M and w: 

i) The relationship between the end-point values m* of 
-^Bondi against M*; see Fig. 10. For M* > 4, we again find 
m* has the same universal value, ~ .864Afpi. 

ii) The relationship of vs a* (once GDH becomes time- 
like). As before, by an appropriate shift, we find a y^j^ that 
can be used as a universal time coordinate for all cases. 

iii) The dependence of F* and M^^^^^ on a* and see 
Fig. 11. For a fixed value of Madm the plots are indis- 
tinguishable, so that even for this broader class of matter 
profiles, there are two universal curves, one for the dynamics 
of and the other for M^^^^^^. In particular, for a given 
w > 0, the time evolution F* and M^^^^^ is identical to that 
obtained with the shell collapse {w = 0). 

In the classical theory, if the collapsing matter is com- 
pactly supported on , to the future of this support the ge- 
ometry is universal, determined by the ADM mass Madm- 
This is because stationary, classical, CGHS black holes are 
characterized completely by A/adm- Whether the situation 
would have a direct counterpart in the semi-classical theory 
is not a priori clear because the semi-classical solutions are 
not stationary and there is no reason to expect the solution to 
be characterized just by one or two parameters to the future 
of the support of f^°\ Our results provide a precise sense in 
which universality does persist. As long as the black hole is 
initially macroscopic and the collapse is prompt, we have : i) 
a universal asymptotic time translation d/dy'y^ (Fig 12); and, 
soon after the formation of the GDH, ii) universal dynamics of 
physical observables with respect to the physical asymptotic 
time y-^. 




Figure 12. Plot of dy~ /Az against the separation in z~ from the 
singularity for various values of M and w with a fixed ADM mass 
M* — 6. The functional dependence y^{z^) determines the 
physics of the outgoing quantum state completely [7, 8]. Coinci- 
dence of these curves in the mean-field theory suggests that the out- 
going quantum state is likely to be universal within the class of initial 
data with the same ADM mass, so long as the collapse is prompt. 



IV. DISCUSSION 

The CGHS model provides a useful arena to explore the 
formation and quantum evaporation of black holes. For, the 
classical action is closely related to that governing the spher- 
ically symmetric gravitational collapse in 4 dimensions and, 
at the same time, the decoupling of matter and dilaton fields 
in the model introduces significant technical simplification. 
However, in this paper, we were not concerned with the full 
quantum theory of the CGHS model. Rather, we restricted 
ourselves to the mean-field equations of [7, 8] and explored 
their implications using high precision numerics. 

A. Viewpoint 

Our analysis was carried out in the same spirit that drove the 
investigation of critical phenomena in classical general rela- 
tivity [5, 6]. There, one takes equations of general relativity 
seriously and shows, for example, that black holes can form 
with arbitrarily small mass. From a more complete physical 
perspective, these black holes would have enormous Hawk- 
ing temperature, whence quantum effects would be crucial. 
To know whether black holes with arbitrarily small masses 
can form in Nature, one cannot really rely on the classical 
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Einstein equations. The viewpoint in those investigations was 
rather that, since general relativity is a self-contained, well de- 
fined theory, it is interesting to explore what it has to say about 
such conceptual issues. The results of those explorations led 
to the discovery of critical behavior in gravitational collapse, 
which is of great interest from a theoretical and mathematical 
physics perspective. In the same vein, in the CGHS model, it 
is conceivable [8] that the relative quantum fluctuations of op- 
erators 9, $, may become of order 1 once the horizon mass 
is of the order of, say, ^/M*M-p\. ' ' Suppose this were to 
happen at a point p on the GDH. Then, to the future of the 
null ray from p to T^, solutions 0, $ to the mean-field equa- 
tions discussed in this paper would be poor approximations of 
the expectation values of 0, $ that result from full quantum 
equations. That is, our solutions to the mean-field equations 
would not be physically reliable in this future region. The 
scope of this paper did not include this issue of the physi- 
cal domain of validity of the mean-field approximation. As in 
much of the literature on semi-classical gravity, we considered 
the entire space-time domain where the mean-field equations 
have unambiguous solutions. And as in investigations of criti- 
cal phenomena, our focus was on exploring non-trivial conse- 
quences of these equations. Specifically, we wished to explore 
two questions; Are standard expectations about predictions of 
semi- classical gravity borne out? Do the mean-field dynamics 
exhibit any universal features? 

B. Summary and conjectures in geometric analysis 

We found that some of the standard expectations of semi- 
classical gravity are indeed borne out: The semi-classical 
space-time is asymptotically flat at as in the classical the- 
ory, but in contrast to the classical case is now incomplete. 
Thus, the expectation [25] that the full quantum space-time 
would be an extension of the semi-classical one is viable. 

However, a number of other expectations underlying the 
standard evaporation paradigm turned out to be incorrect. 
Specifically: 

a) The traditional Bondi mass A/g^fj^j is large and negative at 
the end of the semi-classical evaporation rather than of Planck 
size and positive; 

b) The recently introduced Bondi mass ^/b^"^; remains posi- 
tive but is large, rather than of Planck size at the end of evap- 
oration; 

c) The energy flux F^^^ of quantum radiation deviates from 
the Hawking flux corresponding to thermal radiation even 
when the black hole is macroscopic, the deviation becoming 
larger as the evaporation progresses; and, 

d) Along the 'last ray' from the end of the singularity to 1^ , 
curvature remains finite; there is no 'thunderbolt singularity' 
in the metric extending to . 



Note incidentally that in 4-dimensions, when a black hole with Madm = 
Mq has shrunk down through quantum radiation to mass ^A/adm-'^^pi. 
its horizon radius is less than a fermi, and for a super-massive black hole 
with Madm = 10^ Mq, this radius is a tenth of an angstrom. 



The analysis also brought out some unforeseen universali- 
ties. The most striking among them are: 

i) If M* = Mabm/N is macroscopic, at the end of semi- 
classical evaporation jti* -^^BcJndi/-^ assumes a universal 
value, m* w .864A/pi; 

ii) As long as M* is greater than Mp\, there is a universal 
relation: m* = a(l - e-^(*^*)"' ) Mpi, with a « 0.864, 
/? « 1.42, 7 « 1.15; 

iii) An appropriately defined affine parameter y~y^ along is 
a universal function of the area ogdh of the generalized dy- 
namical horizon; 

iv) The evolution of the Bondi mass Mg^^; with respect to 
an invariantly defined time parameter acDH (or y~^^) follows a 
universal curve (and same is true for the energy flux F^^^). 

These results bring out a point that has not drawn the atten- 
tion it deserves: the number of fields plays an important 
role in distinguishing between macroscopic and Planck size 
quantities. If semi-classical gravity is to be valid in an in- 
teresting regime, we must have ^ 1 and the ADM mass 
and horizon area are macroscopic if Madm/N > AGhAIpi 
and a/iV > Gh. (By contrast, it has generally been assumed 
that the external field approximation should hold so long as 
^^ADM > ^^Pi or a > Gh.) Of course the ADM masses can 
be much larger and for analogs of astrophysical black holes 
we would have il/ADM/(-^^^pi) ^ Gh. After a brief tran- 
sient period around the time the GDH is born, dynamics of 
various physical quantities exhibit universal behavior till the 
horizon area a goes to zero. If 7\/adm/(^-^^pi) ^ 1, the uni- 
versal behavior spans a huge interval of time, as measured by 
the physical affine parameter y^^_^ on or the horizon area a. 

All these features are direct consequences of the dynamical 
equations (2.1) and (2.2) for infalling profiles (3.8) charac- 
terized by two parameters M, w. Of course, with numerical 
analysis one cannot exhaustively cover the full range of so- 
lutions, and given the complete freedom to specify incoming 
flux from one can always construct initial data that will 
not exhibit our universal dynamics — for example, after the 
GDH is formed, send in a steady stream of energy with mag- 
nitude comparable to F^^^ . Here we have restricted atten- 
tion to initial data for which the GDH forms promptly, and is 
then left to decay quantum mechanically without further in- 
tervention. Our intuition is that universality is associated with 
a pure quantum decay of a GDH, pure in the sense that the 
decay is uncontaminated by continued infall of classical mat- 
ter carrying positive energy. Therefore, we conjecture that for 
macroscopic black holes formed by smooth infalling matter 
profiles of compact support, these universalities will continue 
to hold soon after the GDH turns time-like. More generally, 
for profiles in which the positive energy flux carried across the 
GDH by the classical fields is negligible compared to the 
negative quantum flux to the future of some ray z+ = z+, the 
universality should also hold in the future region z+ > z+. 

This scenario provides a number of concrete and interesting 
problems for the geometric analysis community. Start with 
initial data (2.6), (2.7) at with /i'' = and a smooth pro- 

(i) 

file /° with compact support for each of the N fields . 
Evolve them using (2.1) and (2.2). Then, we are led to con- 
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jecture that the resulting solution has the following properties: 

1) The maximal solution is asymptotically flat at right future 
null infinity ; 

2) is future incomplete; 

3) A positive mass theorem holds: The Bondi mass M^J^^^ is 
non-negative everywhere on I^; 

4) So long as Madm > N ^/hjG n, the final Bondi 
mass (evaluated at the last ray) is given by AfgJJj^^; « 

5) Fix a solution So with 7\/adm = ^lo ^ No^Jh/G n 
and consider the curve Cq describing the time evolution of 
the Bondi mass in the acDH/^o - M-Qondi/No plane it de- 
fines. Then the corresponding curve c for a solution with 
M /N < Mo/No coincides with Cq soon after its GDH be- 
comes time-like. 



C. Quantum Theory 

Although the mean-field approximation ignores quantum 
fluctuations of geometry, nonetheless our results provide some 
intuition on what is likely to happen near in the full quan- 
tum theory. First, because there is no thunderbolt singularity 
along the last ray, the semi-classical solution admits exten- 
sions in a large neighborhood of to the future of the last 
ray. In the mean-field approximation the extension is ambigu- 
ous because of the presence of a singularity along which the 
metric is C" but not C^. But it is plausible that these ambi- 
guities will be resolved in the full quantum theory and there 
is some evidence supporting this expectation [8, 28]. What 
features would this quantum extension have? Recall that the 
model has N scalar fields and the black hole emits quantum 
radiation in each of these channels. The Bondi mass that is 
left over at the last ray is MBondi ~ 0.8647VMpi. So we have 
(0.864/24)i\/pi units of mass left over per channel. It is gen- 
erally assumed that this tiny remainder will be quickly radi- 
ated away across 1^ , the right future null infinity of the quan- 
tum space-time that extends beyond the last ray. Suppose it is 
radiated in a finite affine time. Then, there is a point p on 1^ 
beyond which Mg^^^j and F^'^^ both vanish. The expres- 
sion (2.16) of F^^^ now implies that is 'as long as' I^. 
This is sufficient to conclude that the vacuum state (of right 
moving fields /i'') on evolves to a pure state on 2"^ (be- 
cause there are no modes to trace over). This is precisely the 
paradigm proposed in [7]. Thus, the semi-classical results ob- 
tained in this paper provide concrete support for that paradigm 
and re-enforces the analogous 4-dimensional paradigm of [29] 
(which was later shown to be borne out also in the asymptoti- 
cally AdS context in string theory [30]). 

All our analysis was restricted to the 2-dimensional, CGHS 
black holes. As we mentioned in section I, while they mimic 
several features of 4-dimensional black holes formed by a 
spherical symmetric collapse of scalar fields, there are also 
some key differences. We will conclude with a list of the 
most important of these differences and briefly discuss their 
consequences. (For a more detailed discussion, see [8].) 

First, in 2-dimensions, surface gravity n and hence, in the 



external field approximation, the Hawking temperature Tnaw, 
is a constant of the theory; it does not depend on the specific 
black hole under consideration. In 4-dimensions, by contrast, 
K and Tnaw depend on the black hole. In the spherical case, 
they go inversely as the mass so one is led to conclude that the 
black hole gets hotter as it evaporates. A second important dif- 
ference is that, in the CGHS black hole, matter fields f'*^ are 
decoupled from the dilaton and their propagation is therefore 
decoupled from the dynamics of the geometric sector. This 
then implies that the right and left moving modes do not talk 
to one another. In 4 dimensions, the /'*^ are directly coupled 
to the dilaton and their dynamics cannot be neatly separated 
from those of geometric fields $, 6. Hence technically the 
problem is much more difficult. Finally, in 4 dimensions there 
is only one 1+ and only one while in 2 dimensions each 
of them has two distinct components, right and left. Concep- 
tually, this difference is extremely important. In 2 dimensions 
the infalling matter is only in the plus modes, f'"_^\ and its ini- 
tial state is specified just on while the outgoing quantum 
radiation refers to the minus modes, /_ , and its final state has 
support only on . In 4-dimension, there is no such clean 
separation. 

What are the implications of these differences? 

Because of the first two differences, analysis of CGHS 
black holes is technically simpler and this simplicity brings 
out some features of the evaporation process that can be 
masked by technical complications in 4 dimensions. For in- 
stance, since the Hawking temperature Tnaw is an absolute 
constant {hn) for CGHS black holes, the standard paradigm 
that the quantum radiation is thermal till the black hole has 
shrunk to Planck size leads to a clean prediction that the en- 
ergy flux should be constant, F^^-w ^ ;jK2^^g tested this 
simple prediction in the mean field approximation and found 
that it does not hold even when the horizon area is macro- 
scopic. In 4 dimensions, since the temperature varies as the 
black hole evaporates, testing the standard paradigm is much 
more delicate. Similarly, thanks to the underlying technical 
simplicity in the CGHS case, we were able to discover scaling 
properties and universalities. We believe that some of them, 
such as the 'end point universality', will have counterparts in 
4 dimensions but they will be harder to unravel. The CGHS 
results provide hints to uncover them. 

The third difference has deeper conceptual implications 
which we will now discuss in some detail. In 4-dimensions, 
since there is a single and a single I+, unitarity of the 
quantum S-matrix immediately implies that all the informa- 
tion in the incoming state can be recovered in the outgoing 
state. In 2 dimensions, on the other hand, there are two dis- 
tinct questions: i) is the S-matrix from XjJ to unitary? and 
ii) is the information about the infalling matter on recov- 
ered in the outgoing state at ? As discussed above, results 
of this paper strongly support the paradigm of [7, 8] in which 
the answer to the first question is in the affirmative; informa- 
tion on is faithfully recovered on I^. However, this does 
not imply that all the infalling information at is imprinted 
on the outgoing state at . 

In the early CGHS literature, this second issue was often 
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mixed with the first one. Because it was assumed that all (or 
at least most) of the ADM mass is evaporated away through 
quantum radiation, it seemed natural to consider seriously the 
possibility that all the information in the infalling matter at 
can be recovered from the outgoing quantum state at X^. The 
key question was then to find mechanisms that make it pos- 
sible to transfer the information in the right-moving infalling 
modes f^'^ to the left-moving modes f^'^ going out to 2^. In 
[14], for example, the 2 dimensional Schwinger model with a 
position dependent coupling constant was discussed in some 
detail to suggest a possible mechanism. 

However, our universality results strongly suggest that 
these attempts were misdirected. The physical content of the 
outgoing quantum state is encoded entirely in the function 
Vsh^^ ) ['''' ^] onip , the right future null infinity of the quan- 
tum extension of the semi-classical space-time. In the family 
of profile functions /|°' we analyzed in detail, the function 
y~Y^{z~) on Ij^ has universal behavior, determined just by the 
total ADM mass. Since only a tiny fraction of Planck mass is 
radiated per channel in the portion of that is not already 
in , it seems highly unlikely that the remaining informa- 
tion can be encoded in the functional form of y~Y^{z~) in that 
portion. Thus, at least for large M* we expect the answer to 
question ii) to be in the negative: information contained in 
the general infalling matter profile on will not be fully re- 
covered at Ij^ . From our perspective, this is not surprising be- 



cause the structure of null infinity in the CGHS model is rather 
peculiar from the standpoint of 4 dimensions where much of 
our intuition is rooted. In 2 dimensional models, is not the 
full future boundary of space-time. Yet discussions of CGHS 
black holes generally ignore I"^ because, as we saw in sec- 
tion 1, even in the classical theory the black hole interpretation 
holds only with reference to Ip . Indeed, for this reason, in- 
vestigations of quantum CGHS black holes have generally fo- 
cused on the Hawking effect and question i) of unitarity, both 
of which involve dynamics only of /_ for which IjJ does 
effectively serve as the complete future boundary. 

In 4 dimensions, the situation is qualitatively different in 
this regard: in particular, the outgoing state is specified on 
all of future null infinity I+, not just on half of it. There- 
fore, if the singularity is resolved in the full quantum theory, 
would be the complete future boundary of the quantum 
space-time and there would be no obstruction for the S matrix 
to be unitary and hence for the full information on to be 
imprinted in the outgoing state on I+. 



Acknowledgements: We would Uke to thank Amos Ori, 
Madhavan Varadarajan for discussions and for their com- 
ments on our manuscript. This work was supported by 
the NSF grants PHY-0745779, PHY-0854743, the Eberly re- 
search funds of Penn State, and the Alfred P. Sloan Founda- 
tion. 



[1] C. G. Callan, S. B. Giddings, J. A. Harvey and A. Strominger, 
Evanescent black holes, Phys. Rev. D 45, R1005-R1009 (1992). 

[2] Quantum mechanics of black holes, S. B. Giddings, 
arXiv : hep-th/ 94 12 13 8; A. Strominger, les Houches 
lectures on black holes, arXiv : hep-th/ 9501071 

[3] G. Mandal, A. M. Sengupta and S. R. Wadia, Classical solu- 
tions of 2-dimensional string theory. Modern Physics Letters 
A6, 1685-1696 (1991) 

E. Witten, String theory and black holes, Phys. Rev. D 44, 314- 
324(1991) 

[4] A. Ashtekar, F. Pretorius and F. M. Ramazanoglu, Sur- 
prises in the Evaporation of 2-Dimensional Black Holes, 
arXiv: 1011 : 6442vl 

[5] M. Choptuik, Universality and scaling In Gravitational collapse 
of a massless scalar field, Phys. Rev. Lett. 70, 9 (1993). 

[6] C. Gundlach and J. M. Martin-Garcia, Critical phenomena in 
gravitational collapse, LivingRev. Rel. 10:5 (2007) 

[7] A. Ashtekar, V. Taveras and M. Varadarajan, Infornmation is 
not lost in the evaporation of 2-dimensional black holes, Phys. 
Rev. Lett. 100, 211302 (2008) 

[8] M. Varadarajan, A. Ashtekar and V. Taveras (pre-print, 2010) 

[9] F. M. Ramazanoglu and F. Pretorius, Two-dimensional quan- 
tum black holes: Numerical methods. Class. Quantum Grav. 
27, 245027 (2010) 
[10] R. P. Geroch and G. T. Horowitz, asymptotically simple does 
not imply asymptotically Minkowskian, Phys. Rev. Lett. 40, 
203-206 (1978) 

[11] A. Ori, Approximate solution to the CGHS field equations for 
two-dimensional evaporating black holes, Phys. Rev. D 82, 



104009 (2010) 

[12] T. Piran and A. Strominger, Numerical analysis of black hole 
evaporation, Phys. Rev. D 48, 4729-4734 (1993) 

[13] H. Bondi, M. G. J. Van den Berger, and A. W. K. Metzner, 
Gravitational waves in general relativity VIII. Waves from axi- 
symmetric islated systems, Proc. Roy. Soc. (London) A269 21- 
52 (1962) 

[14] L. Susskind and L. Thorlacius, Hawking radiation and back re- 
action, Nucl. Phys. B 382, 123-147 (1992) 

[15] S. Hayward, Cosmic censorship in 2-dimensional dilaton grav- 
ity, Class. Quant. Grav. 10, 985-994 (1993). 

[16] A. Ashtekar and B. Krishnan, Isolated and dynamical horizons 
and their applications, LivingRev. Rel. 7:10 (2004) 

[17] A. Ashtekar and B. Krishnan, Dynamical horizons: Energy, an- 
gular momentum and balance laws, Phys. Rev. Lett. 89, 261 101 
(2002); 

Dynamical horizons and their properties, Phys. Rev. D 68 
104030 (2003) 

[18] D. A. Lowe, Semiclassical approach to black hole evaporation, 
Phys. Rev. D 47, 2446-2453 (1993) 

[19] T. Tada and S. Uehara, Consequences of Hawking radiation 
from 2-d dilaton black holes, Phys. Rev. D 51 4259-4264 ( 1 995) 

[20] S. W. Hawking and J. M. Stewart, Naked and thunderbolt singu- 
larities in black hole evaporation, Nucl. Phys. B 400, 393-415 
(1993) 

[21] D. Yeom and H. Zoe, Semi-classical black holes 
with large N re-scaling and information loss problem 
arXiv: 0907 . 0677vl . 

[22] J. G. Russo, L. Susskind and L. Thorlacius, The Endpoint of 



18 



Hawking radiation, Phys. Rev. D 46, 3444 (1992) 

[23] A. Bilal and C. G. . Callan, Liouville models of black hole evap- 
oration, Nucl. Phys. B 394, 73 (1993) 

[24] A. Strominger, Unitary rules for black hole evaporation, 
arXiv:hep-th/9410187. 

[25] S. W. Hawking, particle creation by black holes, commun. 
Math. Phys. 43, 199-220 (1975) 

[26] S. B. Giddings and W. M. Nelson, Quantum emission of 2- 
dimensional black holes, Phys. Rev. D 46, 2486-2496 

[27] T. Vachaspati and D. Stojkovic, Quantum Radiation from Quan- 



tum Gravitational Collapse, Phys. Lett. B 663, 107 (2008) 
[28] D. Levanony and A. Ori, Interior design of a two-dimensional 
semiclassical black hole: Quantum transition across the singu- 
larity, Phys. Rev. D 81, 104036 (2010) (1992) 
[29] A. Ashtekar and M. Bojowald, Black hole evaporation: A 

paradigm. Class. Quant. Grav. 22, 3349-3362 (2005) 
[30] G. Horowitz, A. Lawrence and E. Silverstein, Insightful D- 
branes, JHEP 0907:057 (2009) 



